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Abstract 

The  need  to  validate  and  improve  numerical 
models  of  the  MPD  thruster  flow,  using  detailed 
comparisons  with  experimentally  determined  plasma 
properties,  is  essential  in  the  effort  to  improve  MPD 
thruster  efficiency.  This  research  represents  an  attempt  to 
experimentally  identify  trends,  in  the  variation  of  plasma 
properties  inside  of  a  self-field  MPD  thruster,  which  could 
be  meaningfully  compared  with,  numerical  results. 
Experimental  results,  consistent  with  the  numerical  model 
in  trend,  show  an  axial  variation  of  the  electron 
temperature  along  the  anode  which  exhibits  a  very 
distinctive  double-peaked  profile.  Experimental  axial 
profiles  of  the  electron  number  density  and  current  density 
were  also  consistent  with  the  numerical  model  in  trend. 
Such  evidence,  along  with  results  of  the  numerical  model, 
suggest  that  an  oblique  shock  is  present  near  the  inlet,  and 
is  required  to  establish  the  radial  pressure  gradient  so 
prevalent  in  the  MPD  thruster.  This  work  marks  one  of 
the  first  non-trivial,  successful  comparisons  between  a 
numerical  model  and  experimentally  determined  plasma 
properties  inside  the  MPD  thruster. 

I.  Introduction 

The  magnetoplasmadynamic  (MPD)  thruster  is  a 
space  propulsion  device  with  the  potential  to  significantly 
enhance  the  payload  mass  delivered  to  GEO  and/or  greatly 
reduce  launch  costs  by  enabling  launch  vehicle 
downsizing.  A  recent  study  has  shown  that  the  MPD 
thruster,  due  to  its  higher  specific  impulse,  is  capable  of 
providing  the  above  benefits  in  excess  of  those  provided 
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by  the  hydrogen  arcjet,  with  only  a  small  increase  in 
transfer  time[l].  The  MPD  thruster  operates  by  using  the 
Lorentz  force  to  accelerate  a  plasma  to  exhaust  velocities 
typically  in  the  10-50  km/sec  range.  Unfortunately,  two 
issues  must  be  resolved  to  actualize  such  benefits: 
thruster  efficiency  and  lifetime.  The  MPD  thruster’s 
insistence  on  operating  inefficiently  has  prohibited  its 
development  for  use  on  spacecraft,  and  has  even  stunted 
current  MPD  thruster  R&D  efforts  in  favor  of  more 
promising  plasma  thrusters.  The  thruster  efficiency  issue 
is  tire  focus  of  this  work. 

Faced  with  the  problem  of  low'  thruster 
efficiency,  MPD  thruster  research  activities  have 
emphasized  the  investigation  of  various  loss  mechanisms 
in  the  thruster[2,3].  As  part  of  this  focus  on  loss 
mechanisms,  there  has  been  a  substantial  increase  in  the 
develpment  of  numerical  codes  to  investigate  the  flow 
physics  inside  of  a  self-field  thruster.  As  discussed  in  a 
1991  review,  only  recently  has  it  been  conceivable  to 
directly  compare  numerical  models  with  experiment^].  - 
initially,  comparisons  of  numerically  calculated 
thrust  versus  terminal  current  curves  compared  quite  well 
to  experimentally  determined  curves  for  self-field  MPD 
thrusters [4,5].  This  success  came  as  no  surprise,  at  least 
in  the  electromagnetic  regime  of  self-field  MPD  thruster 
operation,  because  the  electromagnetic  thrust  is  essentially 
independent  of  propellant  species  and  the  current  density 
distribution  inside  of  the  thruster,  as  noted  in  the 
derivation  of  the  Maecker  thrust  law [6].  Model 

comparisons  With  experimentally  determined  current 
contours  inside  the  thruster  have  also  been  fairly 
successful  qualitatively  and  somewhat  successful 
quantitatively[4,7-ll].  In  the  self-field  MPD  thruster, 
these  comparisons  are  more  relevant  to  the  specific 
impulse  of  the  thruster,  and  do  not  directly  address  the 
efficiency  problem.  It  is  interesting  to  note  the 
experimental  observations  of  Turchi,  et  al.  [12],  who 
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noticed  that  the  current  contours  inside  the  MPD  thruster 
can  be  drastically  altered  by  simply  changing  the  flow 
injection  scheme,  and  that  such  a  change  had  absolutely 
no  effect  on  the  terminal  voltage,  and  most  likely,  no 
effect  on  the  thruster  efficiency. 

The  prediction  of  the  terminal  voltage  is  more 
relevant  to  the  efficiency  problem,  because  the 
identification  of  any  method  to  reduce  the  thruster  voltage 
(for  a  given  current  and  mass  flow  rate)  is  equivalent  to 
increasing  the  efficiency  of  a  self-field  MPD  thruster. 
Until  recently,  numerical  models  had  shown  little  progress 
at  predicting  the  terminal  voltage.  Among  all  of  the 
models,  the  predicted  voltages  were  universally  lower  than 
their  experimental  counterparts.  To  investigate  this 
deficiency,  researchers  have  focused  on  such  loss 
mechanisms  as  microturbulence  effectsflO],  viscous  and 
diffusive  effects[13],  and  the  Hall  effectfl  1,14]  on  thruster 
operation.  Such  efforts  have  been  fruitful,  especially  the 
work  investigating  the  Hall  effect  on  the  anode  fall 
voltage[ll,14].  In  that  work,  terminal  voltages  predicted 
by  a  numerical  model,  which  included  the  Hall  effect  on 
the  voltage  drop  near  the  anode,  compared  very  well  with 
experiment  in  both  magnitude  and  trend. 

Comparisons  of  numerical  results  with  terminal 
and  thrust  measurements  are  not  enough  to  validate  and 
improve  numerical  models[3].  Research  efforts  are  now 
shifting  towards  the  detailed  comparison  of  numerical 
results  with  experimentally  determined  plasma  properties 
in  the  MPD  thruster  interior[15,16].  Trivial  comparisons, 
of  trend  and  magnitude,  are  easily  made  with  the  most 
rudimentary  models,  such  as  the  electron  density'  variation 
with  radius,  or  the  peak  in  electron  temperature  off  the 
cathode  tip.  Prior  to  this  investigation,  there  has  been 
only  one  non-trivial,  satisfying  comparison  between 
numerical  model  and  measurements  taken  inside  the  MPD 
thruster  (a  comparison  of  Niewood’s  numerical  modelfll] 
to  floating  voltage  measurements[17]  at  a  single  location 
near  the  anode  at  various  thruster  currents). 

As  models  become  increasingly  more  powerful 
and  predictive,  the  move  to  the  model  validation  phase 
will  probably  be  hindered  by  the  lack  of  high-fidelity 
plasma  property  data  taken  inside  the  MPD  thruster. 
Arguably  the  most  extensive  study  of  internal  plasma 
properties,  that  of  Tahara,  et  al.[18],  is  of  limited  use  in 
the  further  refinement  of  present-day  models  because  most 
of  the  experimentally  observed  trends  are  dominated  by  the 
error  in  the  experimental  techniques.  These  techniques, 
almost  exclusively  Langmuir  probes,  magnetic  field 
probes  and  emission  spectroscopy [3],  have  been  (and  still 
are)  instrumental  in  the  investigation  of  the  plasma 
physics  associated  with'  the  MPD  thruster.  However, 
plasma  property  data,  internal  to  the  MPD  thruster,  that 
can  be  meaningfully  compared  with  numerical  models  are 
practically  non-existent. 

The  objective  of  this  research  was  to  attempt  to 
experimentally  identify  trends,  in  the  variation  of  plasma 
properties  inside  of  the  MPD  thruster,  which  could  be 
meaningfully  compared  with  numerical  results.  In 
particular,  for  this  study,  the  numerical  results  of 
Niewood[l  1]  were  used  for  comparison  to  experimental 
results.  The  approach  was  to  investigate  the  axial 
variation  of  plasma  properties  in  the  anode  region  of  a 


self-field  cylindrical  MPD  thruster.  A  self-field  cylindrical 
MPD  thruster  was  used  because  this  configuration  is 
conducive  to  comparisons  with  numerical  models.  Probe 
diagnostics  were  utilized  because  of  their  ease  of  use,  and 
the  presumption  that  the  systematic  error  associated  with 
the  Langmuir  probe  wms  acceptable  in  the  anode  region. 
Finally,  the  anode  region  was  chosen  because  it  is  very 
relevant  to  the  thruster  efficiency  problem[ll,14,19],  and 
because  this  region  was  the  focus  of  the  modeling  effort  of 
Niewood[ll].  Fortunately,  and  almost  serendipitously, 
the  above  approach  was  successful  in  identifying  clear 
trends  in  the  axial  variation  of  plasma  properties  along  the 
anode  that  were  consistent  with  the  numerical  model. 

II.  Experimental  Apparatus 

Experiments  were  performed  with  a  cylindrical 
self-field  MPD  thruster  operated  in  a  quasi-steady  mode  (1- 
2  msec  pulse  duration)  at  the  Electric  Propulsion 
Laboratory  at  the  Phillips  Laboratory  at  Edwards  AFB, 
CA.  A  more  complete  description  of  the  facility,  MPD 
thruster,  and  diagnostics  used  in  this  research  can  be  found 
in  reference  15,  20,  and  21. 

Facility 

The  vacuum  chamber,  2.4  m  in  diameter  and  3.6 
m  in  length,  was  pumped  by  two  25  cm  (10  inch) 
diffusion  pumps  to  vacuum  pressures  in  the  low  10‘4  torr 
range,  which  w’as  the  vacuum  pressure  prior  to  firing  the 
thruster.  Once  ignited,  typically  by  a  1  kV  pulse,  the 
thruster  was  powered  by  a  20  kJ,  800  V  maximum  pulsed 
forming  network  (PEN). 

To  ensure  sufficiently  low  vacuum  during 
thruster  operation,  the  propellant  to  the  thruster  was 
supplied  in  a  pulsed  manner  (typical  pulse  length —60-80 
msec).  A  piezoelectric  pressure  transducer  was  used  to 
monitor  tire  gas  pulse,  and  to  verify  that  the  thruster  was 
ignited  during  steady  flow.  The  uncertainty  in  -the  mass 
flow  rate,  which  w’as  calibrated  in  situ  before  each  series 
of  tests,  was  less  than  2%  (due  primarily  to  variations  in 
the  plenum  pressure).  Research  grade  argon  was  the 
propellant  in  all  testing 

To  investigate  the  effect  of  the  residual  vacuum 
on  thruster  operation,  the  voltage  to  the  thruster  was 
measured,  at  a  fixed  PFN  charging  voltage,  as  the  timing 
between  the  gas  pulse  initiation  and  the  ignition  of  the 
thruster  was  varied.  The  resulting  data  showed  that  when 
the  time  delay  was  too  short,  the  voltage  increased,  and 
voltage  oscillations  were  observed  (characteristic  of 
thruster  operation  at  insufficient  mass  flow).  When  the 
time  delay  was  too  long,  the  thruster  voltage  gradually 
decreased,  as  the  entrainment  of  residual  gas  began  to 
affect  thruster  operation.  In  these  experiments,  the  lime 
delay  was  fixed  in  between  the  above  extremes,  where  the 
terminal  voltage  was  observed  to  be  independent  of  the 
time  delay. 

MPD  Thruster 

The  cylindrical  MPD  thruster  used  in  these  tests 
was  identical  to  that  used  previously  at  the  Phillips 
Laboratory[20,21],  The  copper  anode  was  10.16  cm  OD, 
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7.62  cm  ID  and  4.20  cm  long  (from  the  boron  nitride 
injection  plate),  and-  Lite  thoriated  tungsten  cathode  was 
1.27  cm  in  diameter  and  2.3  cm  long.  Propellant  was 
injected  at  the  back  of  the  thrust  chamber,  through  a  boron 
nitride  injection  plate,  via  a  series  of  16  -  0.32  cm 
diameter  injection  holes  (8  at  a  radius  of  1.6  cm  and  8  at 
2.5  cm)  and  through  an  annulus  with  an  outer  radius  of 
0.95  cm  and  inner  radius  of  0.64  cm  at  the  base  of  the 
cathode. 

Shown  in  Fig.  1  is  the  voltage-current 
characteristic  of  the  thruster  at  0.5  g/sec  of  argon.  All 
measurements  reported  in  this  paper  were  obtained  at  a 
mass  flow  rate  of  0.5  g/sec  and  at  currents  below  the  onset 
current  of  5.6  kA.  The  voltage  of  the  thruster  was 
measured  with  the  help  of  a  Tektronix  1000:1  voltage 
probe  (±5%),  and  the  current  with  a  Pearson  Electronics 
current  monitor  (Model  No.  301X,  3%/msec  droop).  The 
uncertainties  in  the  voltage  and  current  were  dominated  by 
shot-to-shot  variation  (at  low  currents)  and  either  high- 
frequency  fluctuations  or  error  introduced  because  the 
waveform  was  not  quite  a  flat-top  profile  (at  high 
currents).  As  seen  previously  in  this  thruster[20,21], 
onset  occurs  (based  on  the  10%  peak-to-peak  voltage 
fluctuation  criterion)  at  a  relatively  low  current  because 
most  of  the  propellant  is  injected  through  the  cathode 
annulus. 

Diagnostics 

The  measurements  reported  in  this  paper  are  only 
a  small  subset  of  a  large  data  base  of  internal 
measurements  taken  in  this  MPD  thruster.  An  induction 
probe  was  used  to  measure  the  magnetic  field,  and  via 
Maxwell's  equations,  derive  the  current  density.  A  triple 
Langmuir  probe[22]  was'  used  to  obtain  the  electron 
temperature,  Te ,  and  electron  number  density,  ne.  A 
floating  probe,  along  with  the  triple  probe  measurements, 
was  used  to  measure  the  plasma  potential  and  anode  fall 
voltage.  No  plasma  potential  or  anode  fall  data  will  be 
reported  in  this  paper;  consult  reference  15  for  further 
details  on  these  and  other  measurements. 

The  induction  probe  consisted  of  a  75  turn  coil 
(1.6  mm  in  diameter,  1.6  mm  long)  placed  in  a  4  mm  OD 
quartz  tube[15].  The  coil  was  oriented  to  measure  the 
azimuthal  magnetic  field,  which  is  the  only  component  in 
a  self-field  MPD  thruster.  Since  the  induction  probe  is 
responsive  to  dB/dt  only,  an  integrator  was  employed  to 
obtain  the  time  history  of  the  magnetic  field. 

The  induction  probe/integrator  system  was 
calibrated  by  first  placing  the  coil  at  the  largest  radius  and 
as  far  into  the  thrust  chamber  as  possible.  Then  by  firing 
the  thruster,  the  magnetic  field  is  calculated  by  Ampere's 
law  (assuming  that  the  total  thruster  current  is  enclosed  by 
the  coil,  B=pLoJ/2rar),  and  correlated  with  the  voltage  from 
the  induction  probe/integrator  system.  After  firing  the 
thruster  at  many  different  currents,  a  least-squares  fit  is 
used  to  determine  the  proportionality  factor  between  the 
voltage  and  magnetic  field.  The  uncertainty  in  the 
magnetic  field  measurements  was  dominated  by  shot-to- 
shot  variation  (-10-20%). 

A  spatial  map  of  die  magnetic  field  inside  die 
MPD  thruster  was  obtained  at  two  operaung  conditions 


(4.4  and  4.8  kA).  At  each  location,  the  magnetic  field 
was  taken  to  be  the  average  of  the  measurements 
corresponding  to  three  firings  of  the  thruster.  From  the 
magnetic  field  map,  and  Ampere’s  law  (Jencl.=2-7trB/jJ.0)> 
an  enclosed-current  contour  map  is  constructed.  The 
contour  map  of  die  enclosed  current  is  then  fitted  with  a 
3rd  order  polynomial  surface  fit;  this  surface  fit  is  used  to 
interpret  the  triple  probe  measurements. 

The  triple  probe  consisted  of  three  cylindrical 
electrodes  (radius=0.063  mm,  lengtb=5.2  mm)  spaced  1.7 
mm  apart,  aligned  with  the  thruster  axis,  and  configured 
such  that  all  electrodes  were  approximately  at  the  same 
radius'  from  the  thruster  centerline.  With  these 
dimensions,  and  configuration,  the  following  effects  can 
be  neglected:  collisional  effects,  magnetic  field  effects, 
electrode  end  effects,  sheath  interactions  between  the 
electrodes,  and  gradient  effects.  Throughout  testing,  the 
Langmuir  probe  was  ion  bombardment  cleaned  by 
establishing  a  glow  discharge. between  the  probe  and  the 
thruster  anode.  This  cleaning  was  performed  every  10 
shots,  which  was  chosen  as  a  compromise  between  the 
need  to  eliminate  contamination  effects  (cleaning  is  desired 
after  every  measurement)  and  time  constraints.  The  probe 
voltage  and  current,  used  to  calculate  Te  and  ne,  were 
taken  as  the  average  of  three  shots.  As  with  the  floating 
probe,  triple  probe  measurements  were  taken  at  eight 
current  levels,  ranging  from  2.2  to  5.3  kA,  and  12  axial 
locations  along  the  anode.  In  addition,  radial  profiles  at 
three  axial  locations  were  obtained  inside  of  the  thruster  at 
4.4  and  4.8  kA. 

The  electron  temperature  and  number  density 
were  calculated  with  equations  12  and  18  of  reference  22 
assuming  7)  =Te  and  singly  ionized  argon  ions.  These 
expressions  were  derived  from  Laframboise’s  exact 
calculations  for  the  ion  .saturation  current  to  a 
cylinder[23].  A  unique  feature  of  the  error  estimates  in 
this  paper  is  that  the  values  for  many  of  the  unknown 
parameters,  used  to  bound  the  error,  were  taken  from  the 
calculations  of  Niewood's  model[ll].  The  model  results 
for  the  anode  region  suggest  that  the  value  of  T[fZ{Te  was 
bounded  between  1/3  and  2,  and  the  component  of  the  ion 
flow  perpendicular  to  the  triple  probe  is  negligible 
compared  to  the  ion  thermal  velocity.  This  last  finding, 
which  is  not  valid  in  most  regions  of  the  thruster, 
significantly  reduces  the  error  associated  with  triple  probe 
measurements)!  9,22,24].  Another  source  of  error  is  the 
unknown  ionization  state  of  the  plasma.  •  Since  doubly 
ionized  argon  ions  are  common  in  MPD  thrusters,  the 
uncertainty  in  electron  number  density  due  to  unknown 
ionization  ■  state  is  accounted  for  by  equation  24  of 
reference  22  (the  number  fraction  of  doubly  ionized  argon 
is  assumed  to  be  bounded  by  0  and  0.5[22]).  With  the 
above  bounds,  the  uncertainty  in  electron  number  density 
was  dominated  by  the  uncertainty  of  these  parameters. 
The  random  error-  associated  with  shot-to-shot  variation 
(10-15%)  was  small  compared  to  the  above  systematic 
error.  Alternatively,  the  systematic  error  associated  with 
electron  temperature  (due  primarily  to  the  electron  drift 
perpendicular  to  the  probe)  was  small  compared  to  the 
random  error  (±10%).  Finally,  the  uncertainty  in  axial 
location  was  taken  as  one  half  of  the  probe  length. 
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III.  Experimental  Results 


One  of  the  more  intriguing  results  of  this 
research  is  presented  in  Fig.  2.  Shown  is  the  electron 
temperature,  1.4  mm  from  the  anode,  plotted  as  a  function 
of  axial  position  for  the  MPD  thruster  operating  at  4.4 
kA.  Accounting  for  measurement  error,  it  is  clear  that  the 
electron  temperature  exhibits  a  double-peaked  axial  profile, 
with  the  upstream  peak  located  at  1.2  cm  and  the 
A  downstream  peak  at  3.4  cm.  This  phenomenon  was 
observed  at  all  thruster  currents  examined  in  this  study,  as 
low  as  2.2  kA  and  as  high  as  5.3  kA[15],  with  the  peaks 
in  the  roughly  the  same  locations.  To  the  author's 
knowledge,  this  is  the  first  observation  of  such 
phenomena  in  the  MPD  thruster. 

To  help  explain  this  phenomenon,  and  to 
investigate  its  implications,  the  numerical  results  of 
Niewood[l  1]  were  consulted.  The  thruster  geometry'  used 
in  this  theoretical  study  is  shown  in  Fig.  3,  and  is 
representative  of  the  Constant  Area  Channel  (CAC) 
thruster  used  in  previous  testing[17,25].  As  the  name 
implies,  the  thruster  has  a  cylindrical  geometry,  is  1 1  cm 
in  length,  and  has  a  cathode  and  anode  radius  of  5.2  cm 
and  7.2  cm  respectively.  In  the  model,  the  CAC  thruster 
is  assumed  to  be  operating  on  argon,  consistent  with  the 
experimental  work  of  this  paper.'  The  three-fluid,  two- 
temperature,  2-D  axisymmetric  numerical  code  is  one  of 
the  most  advanced  developed  for  the  MPD  thruster,  and 
includes  many  physical  processes[ll,14].  Most  notable 
.  of  the  included  physics  is  the  Hall  effect,  finite-rate 
ionization/recombinalion  kinetics,  and  a  model  for  the 
anode  fall  voltage. 

The  axial  profile  of  the  electron  temperature 
predicted  by  the  model,  in  the  CAC  thruster  at  a  distance 
of  1.6  mm  from  the  anode,  is  shown  in  Fig.  4  (Fig.  7-2 
of  reference  11).  Although  the  operating  condition  (39.0 
kA,  4  g/sec  argon)  and  geometry  of  the  CAC  thruster  are 
quite  different  from  that  used  in  this  experimental  study, 
the  magnitude  and  double-peaked'  trend  is  very  similar  to 
the  experimental  data  shown  in  Fig,  2.  Also  note  that  the 
double-peaked  Te  profile  was  observed  in  one  other 
numerical  modeling  effort[26],  as  will  be  discussed  later. 

To  understand  this  behavior,  it  is  helpful  to 
qualitatively  examine  the  electron  energy  equation. 
Assuming  that  ohmic  dissipation,  is  balanced  by 
..  collisional  transfer  between  electrons  and  ions,  the 
following  relation  can  be  derived: 


3  k(Te-Ti)_ 
M; 


where  j  is  the  current  density,  e  is  the  electron  charge,  k  is 
Boltzman's  constant,  and  M[  is  the  ion  mass.  Assuming 
Ti  is  relatively  constant,  this  expression  states  that  Te  is 
larger  in  regions  of  high  current  density  and  low  electron 
density.  Although  this  expression  neglects  such  effects  as 
heat  conduction  and  convection,  pressure  work,  viscous 
dissipation,  and  tire  transfer  of  electron  energy  to  and  from 


the  internal  modes  of  the  ions,  it  is  expected  that  these 
terms  will  be  second  order  in  this  simple  analysis. 

Shown  in  Fig.  5  is  an  experimentally  determined 
contour  plot  of  the  enclosed  current  in  the  MPD  thruster 
operating  at  4.4  kA.  The  characteristic  that  most  of  the 
current  is  emitted  from  the  base  of  the  cathode  and  attaches 
at  the  downstream  end  of  the  anode  has  been  observed  in 
this  thruster  previously[21].  Another  map,  taken  at  4.8 
kA  was  also  very  similar  to  Fig.  5.  Fig.  5  show's  that  the 
current,  density,  in  the  region  of  the  downstream  Te  peak 
(z=3.4  cm),  is  very'  large  and  in  fact  peaks  in  this  region  (j 
is  proportional  to  the  density  of  the  contour  lines).  Eq.  1 
suggests,  due  to  large  ohmic  heating,  that  this  local 
maximum  in  current  density  is  responsible  for  the  peak  in 
Te  .  Such  a  correlation  of  j  with  Te  has  been  observed 
experimentally[18,19]  and  computationally[26,27]  many 
times  at  the  cathode  tip.  and  anode.  The  dowmstream  Te 
peak  in  the  model,  shown  in  Fig.  4,  also  corresponds  to 
the  large  current  density  at  the  downstream  edge  of  the 
anode  (see  Fig.  6-2  of  reference  1 1). 

The  upstream  Te  peak,  for  both  tire  model  and 
this  experimental  work,  is  in  a  region  of  low  current 
density.  Thus,  from  equation  1,  it  is  presumed  that  tire 
electron  density'  is  low  in  this  region.  This  is  indeed  the 
case  for  both  tire  model  and  experiment.  Shown  in  Fig.  6 
are  the  numerical  results  for  the  axial  variation  of  electron 
number  density,  at  a  distance  of  1 .6  mm  from  the  anode, 
as  a  function  of  thruster  current.  Clearly,  at  39  kA,  a 
minimum  in  electron  number  density  is  seen  2  to  4  cm 
downstream  from  the  backplate.  In  Fig.  7,  the  measured 
axial  profile  of  electron  number  density  at  4.4  kA  also 
exhibits  a  valley  at  tire  same  location  of  the  upstream  Te 
peak  in  Fig.  2.  Again,  this  behavior  was  observed  at 
thruster  currents  ranging  from  2.2  to  5.3  kA[15J.  As 
expected,  the  magnitudes  of  the  electron  number  densities 
in  Fig.’s  6  and  7  do  not  compare,  because  of  tire 
significantly  different  mass  flow  rates  in  the  CAC  thruster 
(4  g/sec)  and  tire  MPD  thruster  used  in  these  experiments 
(0.5  g/sec). 

Before  moving  on,  two  points  will  be  made 
about  Fig.  7.  First,  although  tire  random  error  associated 
with  the  electron  number  density  measurements  is  non- 
negligible,  tire  majority  of  the  error  is  systematic  due  to 
unknown  T[  and  species  concentrations.  The  upper  bound 
on  ne  w>as  calculated  assuming  Tj/Z{T e=l/3  and  singly 
ionized  argon  ions  only..  The  lower  bound  was  estimated 
assuming  Ti/Z[Te=2  and  a  number  -fraction  of  doubly 
ionized  argon  ions  of  0.5.  Since  these  parameters  vary 
considerably  along  the  anode  (see  Fig.  4  for  example),  it 
may  not  be  valid  to  assume  that  the  actual  axial  profile  is 
simply  a  shift  in  the  vertical  direction  of  the 
experimentally  determined  profile.  •  Secondly,  it  w-as 
mentioned  earlier  that  a  double-peaked  Te  profile  along  the 
anode  was  observed  numerically  by  Sleziona,  et  al.[26]. 
Upon  examining  their  results,  it  is  believed  that  both 
peaks  are  due  to  high  current  densities,  and  not  from  a 
density  minimum. 

In  Fig.  6  the  numerical  model  indicates  that  the 
density  minimum  is  persistent  over  a  large  range  of 
thruster  currents.  As  noted  above,  this  was  also  tire  case 
experimentally!  In  addition.  Fig.  6  indicates  that  the 
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electron  density  minimum  shifts  downstream  as  the 
thruster  current  is  increased  (for  a  given  mass  flow  rate). 
Again,  there  is  consistency  between  the  model  and 
experiment.  In  Fig.  8,  it  is  seen  that  at  2.6  kA,  the 
density  minimum  is  located  at  1.2+0.26  cm  downstream 
of  the  backplate,  while  at  5.1  kA,  it  has  moved 
downstream  to  z=l. 8+0.26  cm.  As  seen  in  Fig.  7,  at  4.4 
kA  the  density  minimum  is  located  between  these  two 
extremes.  Also  shown  in  Fig.  8,  denoted  by  the  open 
symbols,  are  the  locations  of  the  electron  temperature 
peaks. 

With  such  agreement  between  experiment  and 
numerical  model;  a  logical  question  now  arises  from  these 
results:  why  is  there  a  density  minimum  at  the  anode, 
near  the  inlet  of  the  MPD  thruster?  The  answer  to  this 
question  is  found  by  further  examination  of  the  numerical 
results  of  Niewood.  Shown  in  Fig  9  is  a  contour  plot  of 
the  radial  velocity  inside  of  the  CAC  thruster.  This  plot 
shows  that  large  negative  radial  velocities  occur  at  the 
inlet  of  the  thruster,  and  near  the  region  of  the  upstream 
Te  peak.  The  cause  of  this  inward  flow  is  that  as  soon  as 
the  propellant  enters  the  thrust  chamber,  it  is  accelerated 
towards  the  cathode  by  the  Lorentz  force  in  the  negative  r- 
direction.  Even  though  the  r-component  of  the  Lorentz 
force  is  much  smaller  than  the  axial  component,  it  is 
initially  unbalanced  due  to  the  lack  of  a  radial  pressure 
gradient  at  the  inlet  of  the  thruster.  The  rush  of  plasma 
towards  the  cathode  reduces  the  electron  number  density 
near  the  anode. 

.Shown  in  Fig  10  are  die  streamlines  and  pressure 
contours  in  die  first  3  cm  downstream  of  the  inlet.  This 
figure  shows  that  as  the  flow  approaches  die  cathode,  it  is 
turned  axially  by  what  appears  to  be  an  oblique  shock. 
The  oblique  shock,  which  is  diffused  in  part  by  viscous 
effects  and  possibly  by  numerical  smearing,  can  be  seen  as 
the  pressure  rise  in  lower  right-hand  side  of  Fig.  10. 
Consequently,  the  density  minimum,  and  thus  the 
.  upstream  Te  peak,  is  formed  by  the  initial  accelerauon  of 
the  plasma  towards  the  cathode  at  the  inlet  and  the  contact 
of  the  oblique  shock  to  the  anode.  Only,  an  advanced  2-D 
numerical  code,  such  as  that  of  Niewood,  is  capable  of 
explaining  such  phenomenon. 

In  summary,  near  the  inlet  of  die  thrust  chamber, 
an  oblique  shock  is  required  to  establish  the  radial  pressure 
gradient  inside  the  MPD  thruster.  Considering  die 
differences  in  geometries  of  die  CAC  thruster  and  the 
experimental  MPD  .thruster,  and  the  wide  span  in 
operating  conditions  invesdgated  in  this  study,  it  appears 
that  this  phenomenon  may  be  characteristic  of  all  self-field 
MPD  thnisters. 

IV.  Conclusions 

The  objective  of  this  research  was  to 
experimentally  identify  trends,  in  the  variation  of  plasma 
properties  inside  of  a  self-field  MPD  thruster,  which  could 
be  meaningfully  compared  with  numerical  results.  The 
approach  to  this  objective  was  to  investigate  the  axial 
variation  of  plasma  properties  in  the  anode  region  and 
compare  them  to  the  numerical  results  of  an  advanced 
numerical  model  developed  at  M.I.T.  Experimental 


results  showed  that  (he  axial  variation  of  the  electron 
temperature  along  the  anode  exhibits  a  very  distinctive 
double-peaked  profile.  This  unique  feature  was  observed  at 
all  thruster  operating  conditions  (2.2  -  5.3  kA,  0.5  g/sec) 
examined  in  this  study,  and  was  also  predicted  by  the 
model.  Experimental  axial  profiles  of  the  electron  number 
density  (showing  a  minimum  at  the  location  of  the 
upstream  Te  peak)  and  current  densities  (showing  a 
maximum  at  the  .  location  of  the  downstream  Te  peak)  are 
also  consistent  with  the  numerical  model  in  trend.  As 
observed  in  previous  studies,  the  downstream  peak  is 
produced  by  large  ohmic  heating  from  the  high  current 
density  at  the  dowmstream  edge  of  the  anode.  The 
upstream  Te  peak  is  a  consequence  of  an  electron  density 
minimum  created  by  the  inlet  fluid  dynamics  of  the 
thruster. 

With  increased  confidence  in  the  numerical 
model,  it  was  further  used  to  investigate  the  implications 
of  such  an  upstream  Te  peak.  Based  on  numerical  results, 
it  is  created  in  part  by  the  unbalanced  radial  Lorentz  force 
at  the  inlet  of  the  thruster,  and  the  subsequent  oblique 
shock  which  is  required  to  establish  the  radial  pressure 
gradient  inside  the  MPD  thruster.  Due  to  the  resilience  of 
this  phenomenon  to  thruster  operating  conditions,  tills 
internal  flow  structure  at  the  inlet  may  be  a  characteristic 
of  all  MPD  thrusters. 

Although  the  agreement  between  experiment  and 
theory  in  all  of  the  above  profiles,  is  in  trend  only,  it 
marks  one  of  the  first  non-trivial,  successful  comparisons 
between  numerical  model  and  experimentally  determined 
plasma  properties  inside  of  the  MPD  thruster. 
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Fig.  3:  The  geometry  of  the  MPD  thruster  used  in.  Niewood's[l  1  ]  numerical  model. 


Fig.  4:  Model  predictions  for  the  species  temperatures,  1.6  mm  from  the  anode,  for  the  CAC  thruster 
operating  at  39  kA,  4  g/see  argon. 
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Fig.  5:  The  enclosed  current  contours  in  the  experimental  MPD  thruster  operating  at  4.4  kA  and  0.5 
g/sec  argon. 

_  I  =  23 A  kA 

_ I  =  27.3 

. I  =  37.2 

_ I  =  35.9 

_ I  =  29.0 


o.'oo  '  b.02  '  o.'o4  '  o.bs  '  o.’os  '  o.io  '  0.22  '  0.24 

Fig.  6:  Model  predictions  for  the  electron  number  density,  1.6  mm  from  the  anode,  for  the  CAC 
thruster  operating  at  4  g/sec  argon  and  various  current  levels. 
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Fig.  7:  An  axial  profile  of  the  measured  electron  number  density  near  the  anode. 
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Fig.  8:  Axial  profiles  of  the  measured  electron  number  density,  1.4  mm  from  the  anode,  at  2.3  and 
5.1  kA.  The  open  symbols  correspond  to  the  location  of  the  peaks  in  7g.  (0.5  g/sec  argon) 
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Fig.  9:  Model  predictions  of  the  radial  velocity  contours  inside  the  CAC  thruster  operating  at  31:2  kA, 
4  g/sec  argon. 


Fig.  10:  Streamlines  and  pressure  contours,  within  3  cm  of  the  inlet,  predicted  by  the  model  for  the 
CAG  thruster  operating  at  31 .2  kA*  4  g/sec  argon.  - 


